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PRELIMINARY STUDY OF A POSSIBLE 
AUTOMATIC LANDING SYSTEM 

By Windsor L. Sherman and Sylvia W. Winfrey 
Langley Research Center 

SUMMARY 

Navigation and control laws for a possible automatic landing system have been inves- 
tigated. The system makes use of data from an inertial table and either an airborne or 
ground radar to generate signals that guide the airplane to a landing. All landing maneu- 
vers take place within a zone that extends 6000 m out from the touchdown point, 4000 m on 
each side of the runway center line, and 540 m high. The results show that the system can 
adequately control the airplane on steep, curved decelerating approaches to a landing that 
takes place with small errors from the desired landing point and desired airplane attitude. 
The system studied would interface well with the scanning beam microwave landing system 
(MLS). The use of this system with the MLS makes it possible to incorporate an independ- 
ent landing monitor. 

INTRODUCTION 

The approach and landing phase of flight operations is one of the more hazardous 
portions of aircraft flight with approximately 50 percent of the accidents occurring during 
this phase. Inclement weather not only increases the hazards of approach and landing but 
decreases the efficiency of air transportation by delays resulting from waiting stacks and 
required use of distant alternate airports. A reliable and efficient automatic landing sys- 
tem that reduces pilot workload during the approach and landing and the effects of incle- 
ment weather would increase the safety and efficiency of air transportation. Military 
aircraft operations would derive similar benefits from such a system. 

Automatic landing systems are not a new idea; several have been developed for cur- 
rent transport aircraft, but their use in day-to-day operations has been minimal. These 
automatic landing systems, references 1 and 2, are designed to work with the present 
instrument landing system (ILS). The basic mode of operation is riding the ILS beam to 
a landing. These systems are characterized by long, straight approaches and shallow 



glide slopes. This type of operation subjects the airplane to possible undesirable external 
disturbances over long periods of time when it lacks maneuver room. In addition, ILS 
beams on which this system depends are subject to distortions from weather conditions 
and over -flights by aircraft. The ILS mode of landing an aircraft also imposes an objec- 
tionable noise footprint on communities around the airport. 

A study of recent developments in airborne radar such as the phased array antenna 
and of the scanning beam microwave landing system (MLS), references 3 and 4, indicated 
that the tools were at hand for the development of a total airborne automatic landing sys- 
tem. This system should be capable of executing steep, curved approaches where the air- 
plane would be under positive maneuver control at all times, shortening the time to land 
and reducing the exposure of nearby communities to objectionable noise. Several investi- 
gations of such an automatic landing system have already been made. The favorable 
results incorporated in reference 5 reinforced the concept of the use of airborne radar in 
automatic landing systems. Reference 6 presents methods for navigation in the terminal 
area and reference 7 discusses a feed forward automatic landing system for use with the 
MLS. The system described in reference 7 can use the MLS information to guide the air- 
plane along steep, curved approaches which are defined by a model trajectory generator. 

The automatic -landing-system concept reported herein is a feedback type system 
and the airplane flies a space curve that is dependent on the heading error with respect to 
the runway and its position with respect to the runway center line. The system uses air- 
borne digital computers to process the information obtained from an inertial measuring 
unit (IMU) and either an airborne radar with runway enhancement or a scanning beam MLS 
into commands that control the state of the airplane and effect a landing. Digital comput- 
ers were assumed because this type of computer is more adaptable to the complicated 
logic required by the autoland system. When the airborne radar is used in conjunction 
with the MLS, it is possible to provide an independent landing monitor which is an FAA 
(Federal Aviation Administration) requirement (ref. 8) for automatic landing systems. 

The study reported herein concentrated on the investigation of the control laws and 
logic required to control the airplane along steep, curved decelerating landing trajectories. 
As the airplane characteristics are not explicit in the control laws, these laws should be 
adaptable to different airplanes without major change. The system should be applicable to 
both civil and military aircraft operations. However, maneuver constraints dictated by 
passenger comfort indicated that civil operation would put a greater tax on the system 
capability; civil operations were therefore used as a base of the study. 

The results are presented in the form of plots for ground and altitude tracks and 
tables of initial and touchdown conditions. Because winds, turbulence, and instrument 
noise were not considered in the study, the results are considered preliminary in nature. 
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SYMBOLS 


The International System of Units is used throughout this report. All angles are in 
radians. 

b wing span 

Cp> drag coefficient 

Cp, lift coefficient 

rolling -moment coefficient 
C m pitching-moment coefficient 

C n yawing -moment coefficient 

Cy side -force coefficient 

c mean aerodynamic chord 

E total error 

g acceleration of gravity 

H altitude 


V 2 

Vs 


moment of inertia about y-^-axis 

moment of inertia about y2~axis 
moment of inertia about yg-axis 
mass 


p airplane roll rate about airplane axis 

q airplane pitch rate about airplane axis 


R radar range 
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R 


% 

O 

k rtd 

r 

S 

T 

AT 

T 

t 


u 

V 

V 


V 


T 


v 


w 


z 


j 


a 




radar range vector 

two component range vector in z- coordinates 

two component vectors from reflector to landing gate in z^ coordinates 
airplane yaw rate about airplane axis 
wing area 
thrust 

change in thrust 
turning- rate vector 
time 

airplane speed along y^-axis 
airplane speed 
airplane velocity vector 
airplane velocity vector, magnitude 
airplane speed along yg-axis 
airplane speed along y^-axis 
radar coordinates 

airplane coordinates, referred to principal body axes 
quasi inertial coordinates 
angle of attack 
angle of sideslip 
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r 


flight -path angle 


8 a 

S f 

6 h 

6 r 

6 sp 

h 

"j 

e 

0 a 

°AA 

B AC 

6 AG 

d e 

0 EE 

0 EG 

r 

0 

0 


aileron deflection angle 
elevator deflection angle 
flap deflection angle 
stabilizer deflection angle 
rudder deflection angle 
spoiler deflection angle 

runway coordinates, referred to desired landing point 

inertial coordinates, translate with airplane 

Euler pitch angle 

radar azimuth angle 

radar azimuth angle in z^ coordinates 

command radar azimuth angle in z^ coordinates 

azimuth angle for ground radar 

radar elevation angle 

ground radar elevation angle in Zj coordinates 

elevation angle for ground radar 

servo time constant 

Euler roll angle 

Euler heading angle 



Gains : 


k l 

k 2 

k 3 = 0.8 
= 3.0 or 
k 5 = 5.0 

k 6 

k 7 = 5.8 
k g = 4.4 


variable gain 

variable gain 

4.8 

variable gain 


kg = 0.3 


k 10 = 


k ll = 


k 12 " L1 


k^ 3 = 10.0 


k^ = 1.0 or 1.5 


Subscripts: 


command 


E 


error 


initial 


coordinate index, 


limit 



max 


maximum 


min minimum 

o output 

r runway 

A dot over a variable indicates differentiation with respect to time. 

THE AUTOMATIC LANDING SYSTEM 

A block diagram of the automatic -landing-system concept that was studied with the 
major system components and feedback paths is shown in figure 1. The system uses 
information obtained from an airborne radar or the scanning beam microwave landing sys- 
tem (MLS) and an inertial measuring unit (IMU) to guide the airplane from a point in space 
-in the vicinity of the airport along a curved descending trajectory to a landing. While 
executing this trajectory, the speed of the airplane is changed from approach speed to 
landing speed. The block diagram shows both the airborne and ground radar of the MLS. 
When the MLS is used, the necessary data are transmitted to the airplane where they are 
processed according to the methods of appendix A before entering the system through the 
coordinate transformation block of figure 1. This dual radar capability can be used to 
provide the system with an independent landing monitor (ILM). The ILM would be imple- 
mented by using both airborne and ground radars and calculating the guidance commands 
in parallel computers, using the same algorithm. Comparison of the two sets of guidance 
commands would indicate whether the system is functioning properly. 

The Radar and Inertial Measuring Unit 

The equations used to represent the radar and IMU are given in appendix B. The 
form of the equations was selected so that radar and instrument noise could be added with 
ease at a later date. The radar equations are for the airborne radar but can be used to 
represent the ground radar as well. The airborne radar was selected because of the view 
angle constraint this type of radar imposes on the system. Generally speaking, the field 
of view of an airborne radar is restricted to cone about the center line of the airplane with 
a semiapex angle of 1.13446 rad. This restriction is critical in azimuth. However, recent 
developments in airborne radar have increased the azimuth look angle to about ±| rad. 

This development in look angle eases flight restrictions and makes the use of airborne 
radar more attractive. It was assumed that the airborne radar would be a range optimized 
phased array X band system. The range would be accurate to ±3 m, and the angular accu- 
racy is ±3.927 x 10" 3 rad. The phased array radar is ideal for a landing system because 
failure of some of the antenna elements slightly degrades the antenna performance, and the 


7 



failure is not catastrophic. In addition, a phased array radar can take over such other 
duties as weather surveillance, midcourse navigation, and collision avoidance during other 
portions of the flight. 

Because of ground return, the airborne radar would function much better if some 
form of runway enhancement were used. A semiactive reflector that introduces a small 
frequency shift in the return signal would permit the correct signal to be distinguished 
from the ground return. Such a reflector was assumed in the study and it was located 
900 m down the runway from the touchdown point and 100 m off the center line. 

The inertial measuring for the autoland system measures linear accelerations 
which are integrated to obtain the velocity vector of the airplane. The IMU also measures 
the orientation of the airplane in space and is equipped with rate gyros and angular accel- 
erometers that measure the angular rates and accelerations about the body reference axes. 
Inasmuch as the radar, through the measurement of the range rate vector, can also be used 
to determine the airplane velocity vector, these units can be used to check each other. 

Coordinate Systems and Transformations 

Five coordinate systems were used in the autoland system; these are 


(1) 

inertial coordinates 

*1 

(2) 

quasi inertial coordinates 

Z i 

(3) 

airplane coordinates 

y 3 

(4) 

radar coordinates 

X j 

(5) 

runway coordinates 



The subscript j in the symbol for the coordinate takes on the values 1, 2, or 3, 
the values denoting a specific axis. When j = 1, the axis corresponds to the x-axis of 
a usual x,y,z coordinate system. Similarly, j = 2 corresponds to the y-axis and 
j = 3, the z-axis. The transformations between these coordinate systems are given in 
appendix C. 

The inertial coordinates do not rotate and are forced to translate with the airplane. 
The 77 ^ -axis has the same direction as the runway center line. The inertial coordinates 
and their relationship to the quasi inertial coordinates and airplane coordinates are shown 
in figure 2(a). The quasi inertial coordinates are rotated by the angle 4/, obtained by a 
positive rotation about the 773 -axis. The Zj-axis now points in the direction of flight. The 
airplane coordinates are related to the quasi inertial system through the two additional 
angles 6, obtained by a positive rotation about the z 2 -axis, and <$, obtained by a positive 
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rotation about the y^-axis. The rotations through the angles \p, 9, and $ are the 
standard Euler angles and are calculated by the formulas in appendix B. 

The radar coordinate system is related to the airplane coordinates by the two 
angles 9 & and 0 . The angle 0 a is obtained by a positive rotation about the y^-axis 
and 0 , by a positive rotation about the y^-axis. The rotations are taken in the order 
given. The angles 0 a and 0 g are calculated by the formulas given in appendix B. The 
range is measured in radar coordinates and is a single component vector directed along 
the Xj-axis and is positive from the airplane to the target. This places the origin for the 
radar coordinates at the center of the radar antenna. The radar coordinates translate 
with the airplane. 

The runway coordinates use the touchdown point as a fixed origin and the is 
positive toward the far end of the runway. These coordinates, except for j = 3, are the 
negative of the inertial system. This system is used only for the reporting of the results. 


The Turn Computer 

The purpose of the turn computer is to provide a bank angle command that is used 
to control the roll of the airplane and the rotation about the y 2 body axis to produce a 
coordinated banked turn that will reduce the heading and runway centering errors to zero. 
These two errors must be brought to zero simultaneously. A typical situation to be han- 
dled by the turn computer is shown in figure 3. The airplane is located off to one side of 
the runway with an incorrect heading. The heading error is \p c - i//, where i(/ c is the 
runway heading and in this study was assumed to be zero. The runway centering error is 
given by 9^q - 0^^. The turn command must reflect the total error to be eliminated 
which is the sum of these individual errors. The total error E may be written as 


E = (*c - *) * k l( e AC - 0 Aa) 


( 1 ) 


where kj is a gain introduced so that the weight given to each part of the total error 
could be varied. The total error was made proportional to a bank angle command <^> c 


<b = 
Y c 


K-^) + k l(' 


0 AC " 0 Aa) 


(2) 


where k 2 is a gain that determines the proportionality between the command bank angle 
and the total error. 


For the condition shown in figure 3, and 0^ are positive, 0^ negative, 
tf/ = * = 0, and the airplane must roll negatively in order to direct the normal acceler- 
ation to turn it in the correct direction to eliminate the error. The heading error for this 

case is negative and the runway centering error is positive; thus, there is a value of kp 
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called k l5 that makes the total error zero. There are two cases to consider: 


< k. 


kj g k^. 

The case 


k l <k l- 


When 


< kj, an error exists when the turn computer is turned 


on at the start of the 


landing maneuvers and |^ c - rjs > k] 


0 AC " 0 AA|‘ 


. As xp - xp is 


(9 AC - 0^^ j is positive, the sign of the total error is correct and equa- 
tion (2) produces a roll command that turns the airplane to reduce the total error to zero. 


Once the initial error has been brought to zero and 


motion of the airplane causes k 1 


) 


c p c returns to zero, the forward 
k l( 0 AC ' 0 AA 


to increase. Because 


This motion 


) 


k l(°AC " 6 

now dominates the total error, the airplane turns away from the runway 
causes xp c - xp to increase and continues until xp c - xp becomes large enough to agai 
dominate the command and turn the airplane back toward the runway. This interaction 
produces a snaking motion that is totally unacceptable and the airplane does not line up 
for a landing because both errors are not satisfied simultaneously. 


The case k^ = k^.- For k^ = k^, the total error is zero and no immediate turn 

takes place. However, the forward motion causes no change in \p - xp but increases 
0 AC " 0 AA anc * ^ kus develops a total error. However, the total error is positive and 
would command a positive roll instead of a negative one. If kg is negative, a roll in the 
proper direction takes place and the subsequent motion satisfies both the heading error 
and the runway centering error and the airplane motion becomes a smooth curve toward 
the target point. When k^ > k^, a similar process takes place. 

The total error that is controlling the airplane has an interesting physical interpre- 
tation. If kj = kj, the total error is zero, no turning takes place, and the curvature is 
zero. However, if k 1 is very large with respect to k^, the total error is large, the 
turn becomes very sharp, and the curvature becomes large. Thus, the total error, which 
represents the amount the tangent of a space curve must be turned to zero the error, is 
related to the instantaneous curvature of a space curve. This relationship between the 
total error and the curvature of a space curve illustrates the critical nature of the gain k^, 
as it is this gain that determines the curvature and the sharpness of the turn. 


Experience with the use of equation (2) showed that better turns could be obtained if 
\p was used as a feedback into the turn command. With this feedback, equation (2) 
becomes 


^c 





( 3 ) 


Equation (3) with kg negative was adopted as the basic guidance law for the algorithm of 
the turn computer. This guidance law can be used by itself or combined with other guid- 
ance laws for calculating turn commands. In an algorithm with combined commands, 
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equation (3) is always used for the final phase of the landing. Both the basic turn algorithm 
and a combined turn algorithm were used in this study. The combined turn algorithm used 
pure heading turns and straight line segments in combination with the basic turn algorithm. 
The details of the basic and combined turn algorithms are given in appendix D. 


The Letdown Computer 


The purpose of the letdown computer is to determine the required flight -path angle 
for the airplane, generate turn coordination information, and provide a pitch rate signal 
for the elevator autopilot. This computer also provides the signal for lowering the flaps. 
The geometry used to determine the flight -path command y £ is shown in figure 4. If 
H > 20, 


T c = -tan 


-1 


H - 20 
R H + R RTD 


and if H < 20, 


(4) 


T'c 


-tan 


H 

275.0 


(5) 


Equation (5) was used to compute y c for flaring the airplane. Because these flight-path 
command angles are continuously calculated and referenced to the present position of the 
airplane, disturbances to the airplane positions are compensated for in the y commands. 
The command y c was used in the following equation to compute a pitch-rate command 
for the elevator autopilot: 

-JoVc-vi® <6) 


1c = - k 5 


' 14 hr 


g tan cp sin cp 

r 0 ) + + 0 ' 


V T 


where y Q is given by equation (F5). This first term in equation (6) is the y error term 
that supplies the principal control input for the letdown. The last term, the integral of the 
y error, is necessary so that the system will have zero steady -state error. The term 


co 


g tan (p c 


(7) 


supplies turn coordination information, and the sin (p Q term converts this information to 
the airplane body axis yg. The commanded bank angle is used instead of the actual bank 
angle to give some anticipation to the turn coordination because the q response of the 
airplane is slow. 

The a term comes from the equation for the turning rate of an airplane velocity 
vector in space, which is 
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(8) 


m 


= _1_ 
V 2 


V x 


dV 

dt 


and has a component q - o' about the yg body axis. This component and co both 
specify the rotational rate about the yg body axis for a coordinated turn. These terms 
may be equated and solved for a pitch rate, which accounts for the a term in equa- 
tion (6). Equation (6) is the principal part of the letdown algorithm used in this study, the 
details of which are given in appendix E. The y c signal from this algorithm was used 
to control the flaps. When y became less than zero or attained a specific value, the 
flaps were lowered. 


Speed Control 

Speed control was achieved through the use of a digital thrust control system. The 
following equation was used to compute change in thrust required to eliminate a speed 
error: 


AT 


k e( v c - v t) - 


where 


( 9 ) 


k 6 = V c (28.131 + 58.9736 f ) (10) 

If the AT calculated by equation (9) is added to the present thrust T, the new total thrust 
is T + AT. The positions in time of T and T + AT were found in table I. The differ- 
ence in the times associated with these two thrusts indicated how long it should take to 
change the thrust from T to T + AT. The change in thrust AT was added linearly 
over this time. After the thrust was changed, the thrust controls were shut down for 4 sec. 
At the end of 4 sec, V-p and u were sampled to determine whether they were within 
the control system deadband. The deadband for was ±1.0 m/sec and for u was 
±0.1 m/sec^. If the errors were within the deadband, no more changes were made. How- 
ever, if one or both errors were outside the deadband, thrust changes were made until the 
error entered the deadbands. 


Autopilots and Servos 

The autopilots shown in the block diagram of figure 1 are, with the exceptions of the 
roll autopilot, self-explanatory. In the roll autopilot, was limited in order to hold 

down the roll rate of the airplane. The servos are also shown in figure 1 and all necessary 
data are given there. The necessary integrations on the digital computer were performed 

using the state variable convolution method. As the typical power servo of a large trans- 
port aircraft has a natural frequency of about 30 Hz and a damping ratio greater than 1.0, a 
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second-order servo with these characteristics can be adequately represented by a first- 
order system. All servos were rate and displacement limited. A rate limit of 
±0.34907 rad/sec was used for all servos. The displacement limit varied with the sur- 
face, the 6 a . limit was ±0.34907 rad, the S e limits were 0.26180 rad and -0.43633 rad, 
the stabilizer limit was ±0.34907 rad, and rudder deflection was limited to ±0.5236 rad. 

The rate and displacement limits were rarely, if ever, encountered. 

The output of the aileron servo 6 aQ did double duty. In addition to actuating the 
ailerons, this variable was also used to actuate the spoilers for roll control (see appen- 
dixes D and F). If 6 aQ exceeded 0.1745 rad, the spoilers were actuated. Because of 
the rapidity of spoiler actuation, no servos were used. 

The Airplane 

The last block of figure 1 represents the airplane. The airplane considered in this 
study was a large, four -engine subsonic jet transport with a mass of 90 719.4 kg. In the 
study, the airplane was represented by the six-degree-of -freedom equations of motion 
referred to principal body axes. Linear aerodynamics were used. The physical and 
aerodynamic characteristics of the airplane are given in appendix F. 

Comments on the Autoland System 

The operation of the autoland system was studied for the initial conditions described 
by the box labeled "autoland system" in figure 5. It was assumed that the airplane was on 
a course that would cause it to intercept the boundaries of this volume and, as the airplane 
entered this volume, the automatic landing system would be turned on and a landing 
effected. It is not necessary to actuate the landing system when the outermost boundary 
is crossed; it may be turned on any time after this boundary is crossed but must be actu- 
ated before the boundary closest to the runway is crossed. The operation of the system 
is not restricted to the volume shown. The only restriction is that the airplane be so 
positioned that good radar data can be obtained. In the case of an airborne radar, this 
would be about 40 km from the touchdown point, and also the reflector must be in the field 
of view of the radar. Several runs were made from greater distances than those shown 
for the autoland system. The same maneuvers as for the normal distances were noted, 
together with long periods of straight flight. 

It is interesting to compare the volume of operation of the autoland system with that 
of the current instrument landing system (ILS) system. This comparison is shown in fig- 
ure 5; while the autoland system uses a rectangular prism for its zone of operations, the 
ILS uses a pencil type volume with very low glide slope. Automatic landing systems have 
been developed for use with the ILS type landing system, and one is described in refer- 
ence 2. These are beam rider type systems and cannot make the steep curved approaches 
of which the autoland system described in this report is capable. 
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RESULTS AND DISCUSSION 


The automatic landing system discussed in the previous section was programed to 
be studied on a digital computer. In the program, all integrations except for the servos 
were done with the Adams -Bashforth prediction method. The state variable convolution 
method was used when performing integration connected with the servos. 

The results of the study of the automatic landing system are summarized in table II. 
The first 26 entries give the touchdown conditions for various initial conditions for an 
initial altitude of 540 m. The remaining entries are cases included to demonstrate changes 
in initial altitude, the use of straight line segments with the turn computer, and other spe- 
cial cases. Figures 6 to 11 illustrate a few of the first 26 cases. Generally speaking, the 
airplane executed two basic types of flight path to the landing, a sweeping curve as shown 
in figure 6 or a compound curve as shown in figure 7. Inasmuch as all trajectories were 
variations of these basic ones, airplane responses are given for these two; all other exam- 
ples are confined to ground tracks as the altitude tracks for the first 26 cases do not show 
much variation except in the glide-slope angle. 

The two basic ground tracks with the altitude tracks and airplane responses are 
shown in figures 6 and 7, which correspond to cases 9 and 17, respectively, of table II. 

The ground tracks are smooth curves that go from the initial position of the airplane to 
the touchdown point which is located at 0,0. Two altitude tracks are shown on figures 6(b) 
and 7(b), the one with the crossmarks being the desired altitude track and the other, the 
actual altitude track of the airplane. The desired altitude track is for reference only and 
does not represent a feed forward command for the airplane. The airplane responses for 
the parameters 0, (3 , a, 9, y, and the heading angle i//, presented in figures 6 and 7 
are smooth and relatively well damped and are considered satisfactory. The lateral accel- 
eration encountered during a typical landing maneuver was well controlled and in most 
cases was equal to or less than O.Olg; however, when the bank angle reached its limit, 
the maximum lateral acceleration was 0.046g. In both cases, the values quoted are the 
peaks of long period oscillations that occur during maneuvers. The control surface 
motions, which are not shown, did not reach their rate or displacement limit; and the 
spoilers, which were coupled to the aileron displacement, rarely actuated. 

Figures 8 to 11 illustrate some interesting conditions that this system can handle. 
Figure 8 shows the ground track for the initial condition ^ = -6000 m, = ”2000 m, 

H = 540 m, and \ (case 1 of table II). For this case, the airplane crossed the cen- 
ter line of the runway before it had completed its turn but was able to reverse itself and 
make a satisfactory landing, showing that the system can handle an overshoot. Figure 9 
is a ground track for the same position, except that the airplane has a heading of - and 
is moving away from the runway center line (case 4 of table II). For this case, the air- 
plane was able to make a satisfactory landing. Based on airborne radar view angles 
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assumed for this system, the system can handle heading angles that permit a course away 
from the runway between - 0,2368 rad and 0.7785 rad; the exact value is, of course, depend 
ent on the initial position of the airplane. Figures 10 and 11 illustrate the ability 

of the system to handle the initial conditions = -3000 m, = ±4000 m, and 

H = 540 m and are cases 24 and 22, respectively, of table II. 

The first 26 cases of table II illustrate the capabilities when the landing system was 
turned on as the airplane crossed the boundaries of the box marked "autoland" in figure 5. 
In the first 20 cases, satisfactory touchdown conditions were obtained. In cases 21 to 26, 
satisfactory touchdown conditions were not always obtained; the satisfactory cases of this 
group, for example, cases 22 and 25, represent limiting cases for the method used for 
determining the gains and k 2 « 

The gain k^ is position dependent and is redetermined each time the airplane 
starts a new turn. For case 9, figure 6, only one turn is required to aline the airplane 
for landing and k^ was calculated only once at t = 0. In case 17, figure 7, two values 
of kj were required. The first value was calculated at the beginning of the landing 
maneuver, t — 0, and the second at about 38 sec. This point is labeled "k^ switch" in 
figures 7(a) and 7(d). Equation (D6) was always used to determine the first value of k^. 
The subsequent values of k^ may be calculated using either equation (D6) or equa- 
tion (Dll). The gain k^ was changed when either the roll angle 4> 0 or heading angle 
rate became zero. Suitable logic based on ip c - and 0^~, - 0^^ prevented gain 
changes after the airplane was alined with the runway and at the start of each turn after 
kj had been set. 

As used in the basic roll algorithm, the gain k 2 is also initial condition dependent 
and, for the most part, satisfactory results were obtained when k 2 was determined by 
using equation (D7). For some of the more extreme initial conditions, cases 21 to 26, 

k 4 

satisfactory landings were not always achieved. Setting the ratio ~ = 2.02 in equa- 

k l 

tion (D7) and applying it to initial conditions for case 23 produced very satisfactory land- 
ing conditions. The roll error at touchdown was little changed, but the heading error 
improved by a factor of 10 and the er ror was improved by a factor of 3. These 
results indicate that the selection of the constant k^ in equation (D7) is important for 
achieving a satisfactory landing. The very limited results available on this problem 
indicate that a change in the value of k^ from about 3.8 to 4.0 when k^ is determined 

for the second time is necessary. Based on the roll -angle time history, the rerun of 
k 4 

case 23 with = 2.023 indicates that this case and cases with similar initial conditions 
k l 

represent limiting conditions for this autoland system and the airplane considered. 
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Optimization methods were not employed to determine k* and k^. It is thought 
that little would be gained by applying such methods to the determination of k^; however, 
optimization methods may offer a fruitful approach to the determination of k^. 

The results for the vertical plane show that in general the airplane was slightly 

pitched up at touchdown and the touchdown point was in error between 39.3 m (short) and 
367 m (long) with respect to the desired touchdown point. Part of this overshoot was 
caused by an initial altitude that was too high for the longitudinal control system that was 
limited to a 0.10422 rad glide slope. When the altitude was reduced, the overshoot was 
dramatically reduced, as can be seen by comparing cases 21 and 27. Two other factors 
also affect the overshoot in this system; the first is the distance along the ground track 
which affects the time to touchdown and the second is the flare computer algorithm. In 
this system, the ground track is determined by turn calculation and varies with the initial 
conditions at the start of the turning and the values calculated for k^ and k^. In some 
cases, the path was such that there was sufficient time and distance to achieve a touch- 
down close to the landing point and in others, not. The overshoots produced by this factor 
are not large, about 50 m, compared with runway length and can be tolerated. 

Equation (5) was used as the algorithm for the flare computer. This algorithm 
affects both the overshoot of the desired touchdown point and the vertical speed at touch- 
down. The vertical speed at touchdown for the landings listed in table II varied from 
0.34 m/sec to 3.2 m/sec. Vertical touchdown speed should not exceed 0.5 to 0.6 m/sec. 


One of the worst cases for overshoot and vertical touchdown speed is case 23. In 
this case the overshoot was 367 m and the vertical touchdown speed was 3.2 m/sec. The 
flare used in this run was controlled by the output of equation (5). In order to improve 
these conditions, a new equation to compute y c during the flare was developed. This 
equation, which is based on predicting the impact by using the flight-path angle at a 20 m 
altitude is 




tan" 


(H tan y 


20 


20 


( 11 ) 


where y 20 is the flight path angle at an altitude of 20 m. Case 23 was rerun using equa- 
tion (11) to control the flare. This run, listed as case 34 of table II, shows great improve- 
ment in the vertical touchdown speed which is now 0.58 m/ sec. However, the overshoot 
was increased by about 80 m, giving a total overshoot of about 450 m. As touchdown 
occurred in the first third of the 3000 m runway, this overshoot from the desired touch- 
down point was not considered critical. 

Cases 28 and 30 were included to show how the system operates at low altitude. The 

initial conditions used were y^ = -3000 m, y 2 = ±4000 m, H = 200 m, and i/a = ±-|. 

For the low altitude condition, the turn control functioned in the usual manner. The 
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longitudinal control maintained constant altitude until y reached a predetermined value 
of -0,073 rad. When this value was reached, the longitudinal control system initiated a 
normal letdown. This low altitude letdown is compared with a normal letdown in figure 12, 
in which altitude versus time traces for cases 21 and 28 are shown. 

As was pointed out in the previous section, the basic turn command can be combined 
with other commands to guide the airplane to a landing. Cases 31 and 32 combine the 
basic roll command with straight-line segments and free turns. The ground track for 
case 31 best illustrates this combination. This ground track, presented in figure 13, 
starts with a turn to a specified heading and is followed by a period of level flight until tj 2 
reaches the value given by equation (DIO). At this point, if/ was set equal to the runway 
heading and k^ was calculated by equation (Dll). Equation (D6) cannot be used with this 
turn algorithm to compute k^. The airplane made a smooth turn onto the runway and 
landed. This combined form of lateral control produced no significant changes in the let- 
down to the runway or in the airplane responses. When compared with case 9, the only 
significant change introduced other than the form of the ground track is a reduction of 
11.4 sec in the time to land. 

The forward speed of the airplane was controlled during the landing maneuver. This 
involved changing from the approach speed of 77.12 m/sec to 72.43 m/sec, the landing 
speed. This was accomplished through the speed control system described in the previous 
section. Typical speed and thrust time histories for the landing maneuvers are shown in 
figure 14(a) for case 27. This case has a time to land of about 80 sec. When longer flight 
times were involved, the only major change was an increase of flight time at the lower 
speed. For the typical case the speed change command and the command to lower the 
flaps from 0.436 rad to 0.872 rad were given at the start of the landing maneuver. A 
small decrease in speed was followed by an increase to about 78.5 m/sec, this surge being 
caused by the airplane pitching over onto the letdown trajectory. The speed surge was 
eliminated in approximately 10 sec and speed was gradually reduced to 73,4 m/sec, the 
top of the deadband, by 64 sec. At 77 sec, the engines were idled for the touchdown. The 
u term in equation (9) is important; without this term, which supplies damping, the speed 
excursions, although smooth, were large, the speed dropping to about 60 m/sec halfway 
through the maneuver followed by an overshoot. The speed and thrust curves for case 28 
of table II are shown in figure 14(b). As in case 27 the command to change the speed and 
flap setting was given at the start of the landing maneuver, t = 0. Because y had to 
be -0.073 rad or less, the letdown did not start until 41 sec had elapsed (see fig. 12(b)). 
During the first 41 sec the speed was reduced to 71.9 m/sec; but as the letdown started, 
a large speed surge occurred that would be intolerable in the landing maneuver. This 
speed variation occurred because of lack of coordination between speed change commands, 
flap lowering commands, and start of letdown. The command sequences were modified so 
that these commands occurred at the same time, and the resulting speed and thrust curves 
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are shown in figure 14(c), case 33 of table II. As can be seen, the unacceptable speed 
variations were eliminated, which illustrated the importance of coordination of speed 
reduction, flaps, and letdown. 

The results presented in this section are all for a bland environment, that is, no 
noise and no winds. Some preliminary results for head winds and cross winds with shear 
have been obtained for a maximum wind speed of 25.80 m/sec. These results show that 
there is no appreciable effect on the system, although a decrab control system would con- 
siderably relieve the touchdown stresses. As pointed out in the introduction, the purpose 
of this study was not to develop an automatic landing system that could be put into an air- 
plane, but to investigate the concept of an airborne autoland system that, when combined 
with a scanning microwave landing system, would provide a self-contained system with an 
independent landing monitor. The system studied provides satisfactory airplane responses 
and landings over a wide range of initial conditions (see table II) and is considered a feasi- 
ble system. Although the system functioned well, the results indicate that more work is 
necessary to refine the methods for determining the gains and kg in the turn algo- 
rithm and the flare algorithm, equation (5) or equation (11), of the letdown subsystem. 
Optimization methods might be of value in such a study. 

CONCLUDING REMARKS 

An automatic landing system for aircraft that takeoff and land in a conventional man- 
ner has been studied. This system obtains its data from an airborne inertial measuring 
unit and a radar that may be either airborne or located on the ground. When the two 
radars are combined, as could be done with the microwave landing system (MLS), a system 
with an independent landing monitor results. This system is capable of executing steep, 
curved approaches while changing speed from approach to landing conditions. Landing 
maneuvers that give satisfactory landings may be started as close as 3000 m from the 
touchdown point and over a wide range of heading angles with respect to the runway. The 
results show that the automatic landing system studied is a feasible system. System gains 
were not optimized for this study; however, the use of optimization methods for gain cal- 
culation and an improved flare algorithm should improve the system response. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 19, 1974. 
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APPENDIX A 


SIGNALS FROM A GROUND LOCATED TRACKING RADAR 


As pointed out in the body of the report, the autoland system can use data from either 
an airborne radar or a tracking radar located on the ground. The model for the system 
discussed in appendix B is based on the use of an airborne radar. Appendix A discusses 
an interface that will permit the use of data from a ground radar. 


It is assumed that the ground radar is referenced to the runway center line so that 
radar azimuth angles can be read 0 to ± 7 r or 0 to 277 rad and converted to the former. To 
insure good data, the airplane should be equipped with a nose reflector so that the return 
will come to a discrete point. The data required from the ground radar are the range, 
range rate, the elevation angle, and the azimuth angle. This data would be transmitted to 
the airplane where it would pass through an inertial smoothing process (see ref. 9). After 
smoothing, the data would have to be preprocessed before going to the airplane navigation 
system. The preprocessing is quite simple; the range R can be used without modifica- 
tion as can the range rate. The azimuth angle for the ground radar 0 ^q is measured 
in a plane that is parallel to the ZpZg plane, and the elevation angle for the ground 
radar 0 EG is measured in a plane perpendicular to this plane. The required 
angles 6 ^ and 0 EE are obtained as follows from 0^ G and 0 E q* 


0 EE " -0 EG 


(Al) 


and 

e AA = e AG - * * A 


(A2) 


The transformations and equations of the guidance system can now be used to obtain the 
information required to guide the airplane to a landing. 
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APPENDIX B 


AIRBORNE RADAR AND INERTIAL MEASURING UNIT 


The radar was modeled by the following equations : 


R = R. + i L dt 
1 ^0 M 


pt 

^ x 2 / \ 

0 a = % + 1 

1 J o 

sec - r - q sin 0 O + p cos 0„ tan 0„ 

c \ ct- ct j c 

pt 1 

\ 

V 

- q cos 0 O + p sin 0 1 dt 

^ R a 


(Bl) 

(B2) 


(B3) 


The range vector is a vector that connects the vehicle with the target, and its direction 
with respect to airplane coordinates is given by 0 and 0 . The range rates R x R x 
and Rjj-g are obtained by transforming the airplane velocity component along and perpen- 
dicular to the range vector. The integrand of equations (B2) and (B3) are 9^ and 0 , 
respectively. 

The angles 0 a and 0 g are the radar angles between the airplane axes and the 
line of sight. The angle 0 a is the azimuth angle and 6 e is the elevation angle. The 
order of rotation is first 0 a and then 0 g . The positive direction is the same as for the 
Euler angles which follow. 


The standard Euler angles which are used in the transformations presented in appen- 
dix C were also used to represent the angular output of an inertial measuring unit. Because 
the angular displacements were not expected to exceed 0.52360 rad, no singularities would 
be experienced in computing the Euler angles. Therefore, the Euler angles could be com- 
puted by integrating their rates. The equations used to compute the Euler angles are 


* = "h + 1 [(« sin (p + r cos cp'j sec 9 
9 = 9^ + ^ ^q cos <fi - r sin <p ) dt 
* = + 1 f + ( q sin <p + r cos </>) tan 


dt 


The integrands of the above equations are xp, 9 , and </>; and i/a 
initial values of these angles. 


(B4) 

(B5) 

(B6) 


and </>. are the 
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APPENDIX C 

COORDINATE TRANSFORMATIONS 

The basic coordinate system used is an inertial system that has the 77^ -axis pointed 
in the same direction as the runway and the 7/3 -axis positive downward. The T^-axis is 
perpendicular to the Vi>V% plane positive in the direction for a right-handed coordinate 
system. This system was fixed in the airplane at the center of mass and translated with 
the airplane. The runway coordinates were centered at the landing point and are 
related to 7^ by 

h = -(”1 - 90 °) 

C 2 = -(» 2 * 10 °) 

^3 = ^3 

and the altitude is H = -£g. The third set of coordinates was the z^ system. This sys- 
tem was rotated from the inertial system by the angle 4/. The transformation between 
the two coordinate systems is 


(Cl) 


The airplane body axis system of coordinates y • was related to the Zj system 
by the following transformation: 


(C 2 ) 


These two transformations, when combined, give the usual transformation between 

inertial and airplane coordinates. The angles i//, 9, and cp of the transformations (Cl) 
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APPENDIX C - Concluded 


and (C2) are the usual Euler angles. The next transformation used is between airplane 

axes and the radar axes. This transformation is 



X 1 

x 2 

x 3 


cos 0„ cos 6_ 

O. G 

-sin 0„ 

<X 

cos 0 O sin 

d, G 


sin 0 a sin d Q 

cos 0„ 

cL 

sin 0 a sin 0 g 


-sin0 e 

0 

cos 0 g 


where 0„ and 0 o are angles that define the center line of the radar beam and are 
measured in the same sense as and 0. The angle 0 a is measured in the Y-^y^ 
plane and 0 g is measured from the YpY 2 plane to the beam center line. The 
angles 0 a and d Q are generated in the radar model. 
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APPENDIX D 


TURN ALGORITHM 


Basic Turn Algorithm 
The basic equation of the turn algorithm is 


^c “ k 2 


Vn ~ & - k; 


P ) + k^6 


AC 



(Dl) 


where i// is the desired or command heading angle, i// is the heading angle of the air- 
plane obtained from equation (B4), and 0 AC and 0 AA were computed from the following 
equations : 


9 AC 


tan 



(D2) 


e 


AA 


tan 



(D3) 


The output of equation (Dl) was limited in order to prevent too large bank angles. The 
limiting was as follows : 


-<t>0 =<Pn = ( t>l 


if , 2 so, 0, = §; if v 2 <o, 0 t = |. 


The important part of the use of equation (Dl) is the calculation of the gains kj and 
k2* Both of these gains are dependent on initial conditions; in fact, kj must be recal- 
culated at the start of turns. To illustrate that point, one value of k^ would be required 
for the case shown in figure 6. For the case shown in figure 7, it is necessary to recalcu- 
late kj at about 40 sec where the airplane starts a turn in the opposite direction. The 
cue for recalculating k^ was either <fi or \p changing sign. The gain k 2 was much 
simpler to set and the method used will be given first. The gain kg was calculated from 
the equation 


k 2 


i7 exp ('K 



(D4) 


with 
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APPENDIX D - Continued 


4.8 if >. and 3.0 if ^ | 


The following method was used to calculate the gain k 


1 ‘ 


R Zl = cos ( K l " &) + t )2 sin ( K 1 " 
r z 2 = V l sin ( Kl -!//) + r) 2 cos (Kj - 


(D5) 


where K 1 = | for i? 2 51 0 and K l = for ^2 < The Rz i and R zg were used 
in equations (D2) and (D3) to compute 0 AC and 0^ which were used to determine k^ 


from the following equation: 


k l = 


~ K 1 - 


6 AC " 0 AA 


(D6) 


This equation was used to calculate kj at the beginning of the landing maneuver ( k^ J 
and when either <fi Q or \f/ changed sign ^ki)„. This determination of kj was subject 

to the following restrictions that are based on experience on the computer: 

(1) If tf/j = ±|, only (k 1 ) 1 is calculated 

(2) If t? 2 >0, (k^ or (k^S 1.483 

(3) If r/ 2 <0, (k-^ or (kj) 2 g 1.524 

(4) If (k 1 ) 2 >(k 1 ) 1 , use kj ^ in place of (kjjg. 

After ( kj ) 2 had been determined for the first time, logic was used to block further com- 
putation of this gain. 

The gain k 2 , which is also initial condition dependent, was calculated by 
k l 

where k 4 = 3.0 if = | or k 4 = 4.8 if |i//.l<£. 


k 2 = r exp ("l^-^l 


(D7) 


Combined Turn Algorithm 

In the combined turn algorithm used in this study, straight line segments were com- 
bined with equation (Dl). At the start of a landing maneuver, equation (Dl) was used with 

ki = 0 and k2 positive and given by 
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APPENDIX D - Continued 



to turn the airplane to a desired heading angle given by 


(D8) 


'P r 


tan 



(D9) 


The exponential in equation (D8) is used to prevent too rapid a buildup of the roll angle and 
the 8.0 is the number of seconds required for the gain to reach its full value. The 8.0 is 
not a fixed number but can be adjusted as required. Logic was incorporated to limit kg 
to 4.0. The N in equation (D9) is the intercept on the runway center line of the new 
course and can, within limits, be set as desired. In this study, N was set at 3000 m. 

The value assigned to N should be less than but should not be much less than 

3000 m; precise limits have not been established. 


After the airplane had turned through the angle \p c , given by equation (D9), the roll 
angle command was set to zero and kept at that value until 


77 2 = 


V + V • i 
. max min ) 

g tan (pj. 


)2 


(D10) 


At this point ip c was set equal to the runway heading and 



- * ~ k 3^ 
9 AC ' 0 AA| 


was calculated by 
(DU) 


where ip and ip are the actual if/ and \p of the airplane calculated by equation (B4) 
and 6 and are calculated from the output of radar algorithm by equations (D2) 

and (D3). In the combined turn algorithm, k^ had to be greater than or equal to 1.2 and 
kg was set equal to 4.0. Equation (Dl) was used to determine <p c . 


Two equations have been presented for the calculation of the gain k^. These are 
equations (D6) and (Dll). Equation (Dll) is the general version and it must be used to 
determine kj for the combined turn algorithm and may be used in the basic turn algo- 
rithm to calculate ^ k^ j . Equation (D6) is a special form of equation (Dll) and was 
necessary in the basic turn algorithm for the calculation of (kj) . 

The output of equation (Dl), 0 was used to compute 


0 E = 



(D12) 


and the condition 


-0.45 = = 0.45 


(D13) 
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APPENDIX D - Concluded 


was imposed. This 0 E was used as indicated in the block diagram (fig, 1) to obtain a 
6 a for the airplane. 


Calculation of Spoiler Deflection Angles 

The spoiler deflection was tied to the aileron deflection through the following 
equations : 

K 6sp<-f> 6 sp= "I 

If 5 a <-0.15708, 5 = -1,08385 - 6.96 a 

<1q bp O 

If -0.15708 s 5 ao i 0.15708, S sp = 0 

If 6 a > 0.15708, 6_„ = 1.08385 - 6.96 a 
«-o S P o 
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APPENDIX E 
LETDOWN ALGORITHM 

The basic command for the pitch autopilot is a pitch rate given by 

«c = - k 5 k 14(>'c- r o) + 

This equation and the individual components were discussed in the main body of the report. 
However, when used in the autoland system, constraints were used to limit certain of the 
terms. These restrictions are: 

If y c > -0.073, y c = 0 

If y c < -0.10472, y c = -0.10472 

If H > 20.0, k 14 = 1.0 

If H s 20.0, k 14 = 1.5 

This command was used as indicated in the block diagram, figure 1, to produce a 5 eQ 
which was fed to the elevator. 

As previously indicated, flaps and speed commands were keyed to the airplane glide 
slope. The most satisfactory method for flap control was 

If y c <0, 5 f = 0.43633 + 0.08727t (E2) 

If 6f S 0.87266, 6 f = 0.87266 

In addition, as indicated in figure 1, the glide slope was used to change speed com- 
mands. When y c became less than zero, a new speed command was given to the system. 
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APPENDIX F 


AIRPLANE DATA 

The airplane used in this study was a typical four -engine jet airliner, with a 
stretched fuselage, currently in use on civil airlines throughout the world. 

Physical Characteristics 


Wing span, b, m 43.30 

Wing area, S, m 2 267.90 

Mean aerodynamic chord, c, m 6.92 

Mass, m, kg 90719,4 

Inertias : 

I Yl , kg-m 2 3.896 X 10 6 

Iy 2 , kg-m 2 9.9333 x 10 6 

I y , kg-m 2 1.3616 x 10 7 


The airplane was powered by four jet engines, each having a maxirnum thrust of 
8.0068 x 10 4 N. 


Aerodynamic Characteristics 
Longitudinal aerodynamics, c.g. at 25 percent M.A.C.- 



8C l 3C 

o 

II 

— — a + 

8 a 96j 

n 

C L 2 , 8C D 


7.03tt 86 f 

8C l 

= 4.87 


5 f - 0.436332) + 0.7 


(b f - 0.436332 j + 0.037 


da 


where 0.7 and 0.037 are the values of C L and Cy, for a = 0.0. 


86 f 

9C D 


0.916732 

0.161574 


86f 

The expressions for C^, C^, 


ac L 

— — , and 

% 


3C n 

— - are for 0.436332 S 6 f § 0.872664. 
96 f 1 
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APPENDIX F - Continued 


ac 


m 


da 

9C 


m 


da 

ac 


m 


d§ 

9C 


m 


36, 


■1.115 


-0.241 sec 


= -0.707 sec 


= -1.017 


ac 


86 


HI = -1.803 


H 


— = -260.68 N-sec-m’ 
9V 


The flaps were accounted for by varying C D and C L as shown in the list of aerody- 
namic parameters. The lift-curve slopes for the two flap positions were almost the same 
therefore, only one value of dCj^/da was used. 

Lateral aerodynamics, c.g. at 25 percent M.A.C.- 


ac y 

967 " 


36 r 
9C 


n 


n 


96 j 
9C 
"96a 
9C Z 

36 a 


9Cy 


9C y 

90 

9Cy 

dip 

B Jh 

3/3 


0.175 

0.114 


-0.105 


-0.00344 

0.0745 


- 0.866 


0.0539 sec 
0.0881 sec 


- 0.210 
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APPENDIX F - Continued 


3C ? 

— — - -0.111 sec 

30 

3C; 

- — - = 0.0614 sec 
30 


— — = 0.173 
3/3 ' 

3C 

— — = -0.0182 sec 
30 

For the spoiler yawing-moment coefficient: 

ff 6 sp £ -0.34907, C ngp - -0.007166 sp 

If -0.34907 < 5 < -0.14835, C n = -0.0123196 SD - 0.0018275 

sp u sn F 


If -0.14835 £ 6 0 = 0.14835, C =0.0 
sp n g p 

If 0.14835 <6 < 0.34907, C n = -0.0123196 SO + 0.0018275 

sp n gp f 

If 5 5 0.34907, C n = -0.007166 SO 

sp ’ n sp 

For the spoiler rolling-moment coefficient: 


If 

If 

If 

If 

If 


5 ^ -0.34907, C ; = -0.025446 sp 

sp ‘-sp r 

-0.34907 <5, <-0.12217, C ; - -0.041076 SD - 0.0049 

sp t g p op 

-0.12217 S 8 ^ 0.12217, C ? = 0.0 

sp 6 sp 

0.12217 < 6 < 0.34907, C / = -0.041076 SD + 0.0049 

sp 6 sp p 

6 S 0.34907, C; = -0.025445 SD 

sp ’ ^ S p bp 


The airplane was simulated by the six-degree-of-freedom equations of rigid body motion. 
The aerodynamic forces and moments were calculated in stability axes and then trans- 
formed to body axes for the computations of airplane motions. At the start of a landing 
run, the airplane was trimmed for straight and level flight, with the flaps set at 
0.43633 rad. 
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APPENDIX F - Concluded 


In addition to the airplane parameters already determined, i//, 6, and <p in appen- 
dix B, it is also necessary to know the angle of attack a, the angle of sideslip /3 , their 
time derivatives, and the flight -path angle y. The following equations were used to com- 
pute these parameters: 


a 


a = 


•“■‘( 5 ) 

(FI) 

uw - wu 
2 2 

(F2) 

Vl + w' 5 



(F3) 

uv - vu 
U 2 + V 2 

(F4) 

definition for y 

is 6 - a; however, this definition is only useful for 


extremely small roll angles because 6 and a are measured in different planes inclined 
from each other by the roll angle. In this study, the definition adopted for y is the angle 
between the velocity vector and a plane parallel to the ground, with y measured in the 
Zj coordinate system; thus, 


y = -tan'^zg/zj (F5) 

and is compatible with the command value of y that is computed by equations (4), (5), 
and (11). 
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TABLE I.- VALUES OF THRUST USED FOR SPEED CONTROL SYSTEM 


Thrust, T 

"f — — — ■ — - 

Time, t 

5.24890151 X 10 3 

0.25 

1.10315896 X 10 4 

.5 

1.69032421 x 10 4 

.75 

2.44652189 x 10 4 

1.0 

3.33616621 X 10 4 

1.25 

4.09236389 x 10 4 

1.4 

4.67063270 x 10 4 

1.5 

5.56027702 X 10 4 

1.625 

7.11715458 x 10 4 

1.75 

9.45247093 x 10 4 

1.875 

1.13429651 x 10 5 

2.0 

1.25217438 x 10 5 

2.125 

1.31222538 x 10 5 

2.25 

1.38784514 x 10 5 

2.5 

1.44567202 x 10 5 

2.75 

1.49015424 x 10 5 

3.0 

1.53241235 x 10 5 

3.25 

1.56577401 x 10 5 

3.5 

1.60047014 x 10 5 

3.75 

1.62360089 x 10 5 

4.0 

1.64584200 X 10 5 

4.25 

1.66808311 X 10 5 

4.5 

1.69032421 x 10 5 

4.75 

1.71212050 x 10 5 

5.0 

1.73391679 x 10 5 

5.25 

1.74592698 x 10 5 

5.5 

1.76816809 x 10 5 

5.75 

1.77928865 x 10 5 

5.9 



Initial condition 


condition 


Case 

? 1 

? 2 

O '. 

1 

H 

t , 

sec 

1 

-6000 

-2000 

17/2 

540 

101.6 

2 



n /4 

540 

90.6 

3 



0 

540 

90.2 

4 



- tr/6 

540 

93.5 

5 


2000 

- 7r/2 

540 

102.6 

6 



- 7r/4 

540 

90.5 

7 



0 

540 

90.2 

8 



n /6 

540 

93.7 

9 


-4000 

n/2 

540 

121.4 

10 



n / 4 . 

540 

110.0 

11 



0 

540 

111.8 

12 


4000 

- n /2 

540 

121.8 

13 



- tt /4 

540 

109.9 

14 



0 

540 

112.4 

15 

-4000 

-4000 

tt /2 

540 

92.6 

16 



it / 4 

540 

88.1 

17 



0 

540 

88.7 

18 


4000 

- n /2 

540 

92.9 

19 



- 17/4 

540 

88.2 

20 



0 

540 

88.8 

21 

-3000 

-4000 

n /2 

540 

80.9 

22 



tt /4 

540 

77.2 

23 



0 

540 

82.6 

24 


4000 

- n /2 

540 

81.1 

25 



- n /4 

540 

77.3 

26 



0 

540 

79.1 

27 


-4000 

17/2 

528 

79.4 

28 



77/2 

200 

81.1 

29 


4000 

- 17/2 

528 

81.8 

30 



- 17/2 

200 

81.6 

b 31 

-6000 

-4000 

n /2 

540 

110.0 

b 32 



0 

540 

104.7 

c 33 

-3000 

-4000 

n /2 

200 

77.8 

d 34 



0 

540 

83,4 


a Second value of is new kj ir 
b Cases 31 and 32 use the combined t 
c Flap and speed change commands j 
d Same as case 23 except flare algor 
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"stabilizer autopilot c Servo 
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6 H t-0.0377' 



\\)UniIL Indicates equipment fixed to airplane 
that senses airplane motions 

(2) l MU, Inertial measuring unit 

(3) a, a, |J and |3 are computed in the l MU 

(4) All servos are rate and displacement limited 


Rudder autopilot 
6 r “ -k 


. Speed control 
AT=k 6 (V c -V!-mu 
k 6 =V c (2S. 131+58. 9735f) 


Servo 
. i 

1 +tD 
* c 0. 03964 


Airplane motions *- 

Figure 1.- Block diagram of the automatic landing system. 
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CO 
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Notes 


1) Runway coordinates are centered at touchdown point ; positive directions 
shown at right 

2) Airplane inertial coordinates are 1 1 to runway coordinates but centered 
in airplane 

3) The z> coordinates are centered in the airplane and rotated away from 
the airplane inertial coordinates by the Euler angle i/i, the heading angle 
of the airplane. The following sketch shows these relationships: 



4) For the runway i]j * i/i = 0 

5) The Vy V 2 plane coincides 
with the Z|,z 2 plane 


Airplane 

location 



Figure 3.- Geometry used for the turn computer. 



O 



Note: Vertical scale and runway width scale (4) times larger than basic scale. All distances in runway 
coordinates 


Figure 5.- Comparison of autoland and current ILS zones of operation. 
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Figure 6.- Airplane and system response for initial condition = -6000 m, ? 2 “ " 4000 m ’ H ~ 540 m ’ ^ 

(case 9 of table II). Tick marks represent 5-sec intervals. 
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(b) Altitude track with respect to the ground. 
Figure 6.- Continued. 
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(c) Heading angle. 
Figure 6.- Continued. 
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(d) Roll angle. 
Figure 6.- Continued 
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(e) Sideslip angle. 
Figure 6.- Continued. 
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(f) Flight -path angle. 
Figure 6.- Continued 
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(g) Angle of attack. 
Figure 6.- Continued. 
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(h) Pitch angle. 
Figure 6.- Concluded 
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(a) Ground track. 

Figure 7.- Airplane and system response for the initial condition ^ = -4000 m, £2 = “4000 m, H = 540 m, 
if/ = 0.0 (case 17 of table II). Tick marks indicate 5-sec intervals. 
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Heading angle, \j/, rad 



(c) Heading angle. 
Figure 7.- Continued. 
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(d) Roll angle. 
Figure 7.- Continued. 
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(e) Sideslip angle. 
Figure 7.- Continued. 
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(f) Flight -path angle. 

Figure 7.- Continued. 
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(g) Angle of attack. 
Figure 7.- Continued 
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(h) Pitch angle. 
Figure 7.- Concluded 




Figure 8.- Ground track for the initial condition ^ = -6000 m, ? 2 - -2000 m, H 540 m, 

of table II). Tick marks indicate 5 -sec intervals. 
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Figure 9.- Ground track for the initial condition = -6000 m, = “2000 m, H = 540 m, if/. = - 

of table II). Tick marks indicate 5-sec intervals. 


cn } =3 



x 10 1 


500 


400 


300 


£ 2 ,m 


200 


100 


Direction of flight 


i 


f 


• jr 


y 


y 




t— L^l---L~l+l 
200 -200 


Desired touchdown point 

i i 

1 -4—1 — t' 




,-K" 


rTT L I 1 L _L J I I I I I I I .1 .1 L J L_.L J — J — I I .1 L_I L_ 

600 -1000 - 1400 -1800 -2200 -2600 -3000 -3400 -3800 


Figure 10.- Ground track for the initial condition = -3000 m, = -4000 m, H = 540 m, ^ = - ~ (case 24 

of table II). Tick marks indicate 5-sec intervals. 
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Figure 11.- Ground track for the initial condition ^ = -3000 m, £ 2 = -4000 m, H = 540 m, | (case 22 

of table II). Tick marks indicate 5-sec intervals. 
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(b) H = 200 m (case 28 of table II). 
Figure 12.- Concluded. 
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Figure 13.- Basic turn algorithm used with a straight line segment ^ = -6000 m, ? 2 - -4000 m, H 540 m, 

= 5 (case 31 of table II). Tick marks indicate 5-sec intervals. 
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(a) Speed changed and flaps lowered at t = 0 (case 27 of table n). 
Figure 14.- Speed and thrust curves. 
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(b) Constant altitude flight for 39 sec. Speed change and flaps lowered at t = 0 (case 28 of table II). 

Figure 14.- Continued. 


OS 

cn 



NASA -Langley, 




r 


Thrust 


O 

c/> 


TD 
CD 
CD 

& :n\ 


60 U 


400 h 


X 


\ 


Speed command 


X^Speed 


POO I 


0 [■[ [ 


I I i l i.J l I L Li 111.1. 1 i U l U UJlM.H.iJ li I LI I i i J l.i !. I I I I L! I [j U I 
10 20 30 40 50 60 

Time, t, sec 



LLU1J 


(c) Constant altitude flight for 39 sec. Flaps lowered with speed change command (case 33 of table II) 

Figure 14.- Concluded. 
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